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Abstract 

This paper introduces Quibbs vl.3, a Java application available for free. (Source 
O code included in the distribution.) Quibbs is a "code generator" for quantum Gibbs 

^ sampling: after the user inputs some files that specify a classical Bayesian network, 

Quibbs outputs a quantum circuit for performing Gibbs sampling of that Bayesian 
O network on a quantum computer. Quibbs implements an algorithm described in 

earlier papers, that combines various apple pie techniques such as: an adaptive fixed- 
point version of Grover's algorithm, Szegedy operators, quantum phase estimation 
and quantum multiplexors. 
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1 Introduction 



We say a unitary operator acting an array of qubits has been compiled if it has been 
expressed as a Sequence of Elementary Operations (SEO), where by elementary oper- 
ations we mean 1 and 2-qubit operations such as CNOTs and single-qubit rotations. 
SEO's are often represented as quantum circuits. 

There exist software, "general quantum compilers" (like Qubiter, discussed 
in Ref. pj), for compiling arbitrary unitary operators (operators that have no a pri- 
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ori known structure). There also exists software, "special purpose quantum compil- 
ers" (like each of the 7 apphcations of QuanSuite, discussed in Refs.[2l |3l H]), for com- 
piling unitary operators that have a very definite, special structure which is known a 
priori. 

This paper introduce^ Quibbs vl.3, a Java application available for free. 
(Source code included in the distribution.) Quibbs is a "code generator" for quan- 
tum Gibbs sampling: after the user inputs some files that specify a classical Bayesian 
network, Quibbs outputs a quantum circuit for performing Gibbs sampling of that 
Bayesian network on a quantum computer. Quibbs is not really a quantum compiler 
(neither general nor special) because, although it generates a quantum circuit like 
the quantum compilers do, it doesn't start with an explicitly stated unitary matrix 
as input. 

Quibbs implements the algorithm of Tucci discussed in Refs. [HI [TJ |8]. The 
quantum circuit generated by Quibbs includes some quantum multiplexors. The Java 
application Multiplexor Expander (see Ref. |5]) allows the user to replace each of those 
multiplexors by a sequence of more elementary gates such as multiply controlled NOTs 
and qubit rotations. Multiplexor Expander is also available for free, including source 
code. 

For an explanation of the mathematical notation used in this paper, see some 
of my previous papers; for instance, Ref. [H] Section 2. 

Throughout this paper, we will often refer to an operator V. V is defined by 
figure 5 of Ref. [6] . We will also use the acronym AFGA (Adaptive Fixed-point Grover 
Algorithm) for the algorithm described in Ref.[8j. 



2 The Control Panel 

Figjl] shows the Control Panel for Quibbs. This is the main and only window of 
Quibbs (except for the occasional error message window). This window is open if and 
only if Quibbs is running. 



^The reason for releasing the first public version of Quibbs with such an odd version number 
is that Quibbs shares many Java classes with other previous Java applications of mine (QuanSuite 
discussed in Refs.[2j [3l |4], QuSAnn discussed in Ref.[5J, and Multiplexor Expander discussed in 
Ref. so I have made the decision to give all these applications a single unified version number. 
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© e e Quibbs 




Ver. 1.3 
Inputs- 



l/O Folder 3nodes 



Do a pre-run ^ Re ad Bayesian Net ^ 



-Staning State- 



A: al(00)0 
B: bl(0)0 
C: cKOlO 



Random StarT ) 



Number Of Probe Bits (for each PE step) 
Number Of Phase Estimation (PEJ Steps 
Maximum Mumber of Crover Steps 
Camma Tolerance (deqs) 
Delta Lambda (degs) 



!_! Omit V Gates (diagnostic) Write Q.Circuit FileT) 
Outputs 



Starting Camma (degs) 1.44fi6Sle+02 

Prob. of Starting State 9.100000e-D2 

Number of Qubits j 
Number of Elem, Ops. 



.//, 



Figure 1: Control Panel of Quibbs 



The Control Panel allows you to enter the following inputs: 

I/O Folder: Enter in this text box the name of a folder. The folder will contain 
Quibbs' input and output files for the particular Bayesian network that you 
are currently considering. The 1/0 folder must be in the same directory as the 
Quibbs application. 

To generate a quantum circuit, the 1/0 folder must contain the following 3 
input files: 

(Inl) parents.txt 
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(In2) states.txt 
(In3) probs.txt 

A detailed description of these 3 input files will be given in the next section. For 
this section, all you need to know is that: The parents . txt file lists the parent 
nodes of each node of the Bayesian net being considered. The states.txt 
file lists the names of the states of each node of the Bayesian net. And the 
probs.txt file gives the probability matrix for each node of the Bayesian net. 
Together, the Inl, In2 and In3 files fully specify the Bayesian network being 
considered. 

In Fig|T| "3nodes" is entered in the I/O Folder text box. A folder called 
"3nodes" comes with the distribution of Quibbs. It contains, among other 
things, Inl, In2, In3 files that specify one possible Bayesian network with 3 
nodes. The Quibbs distribution also comes with 3 other examples of I/O folders. 
These are named "2nodes", "4nodeFullyConnected" and "Asia". 

When you press the Read Bayesian Net button, Quibbs reads files Inl, In2 
and In3. The program then creates data structures that contain complete infor- 
mation about the Bayesian network. Furthermore, Quibbs fills the scrollable list 
in the Starting State grouping with information that specifies "the starting 
state". The starting state is one particular instantiation (i.e., a particular state 
for each node) of the Bayesian network. Each row of the scrollable list names 
a different node, and a particular state of that node. For example, Figjl] shows 
the Quibbs Control Panel immediately after pressing the Read Bayesian 
Net button. In this example, the Bayesian net read in has 3 nodes called A, B 
and C, and the starting state has node A in state al, node B in state hi and 
node C in state cl. 

Suppose y4 is a node of the Bayesian net being considered. And suppose A has 
Na states. Quibbs will give each state of node A a "decimal name" ; that is, a 
number from through A^^ — 1. The "binary name" of a state is the binary 
representation of its decimal name. As shown in FigjT| the scrollable list of the 
Control Panel gives not only the "english name" of the state of each node, but 
also the binary and decimal names of that state. For example, Figjl] informs us 
that state al of node A has binary name (00) and decimal name 0. 

If you press the Random Start button, the starting state inside the scrollable 
list is changed to a randomly generated one. Alternative, you can choose a 
specific state for each node of the Bayesian net by using the Node State 
Menu, the menu immediately to the left of the Random Start button. To 
use the Node State Menu, you select the particular row of the scrollable list 
that you want to change. The Node State Menu mutates to reflect your row 
selection in the scrollable list. You can choose from the menu a particular node 
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state. When you do so, the selected row in the scrollable list changes to reflect 
your menu choice. 

When you press the Do a Pre-run button, Quibbs both reads and writes files: 
it reads files Inl and In2 (but not In3, so if In3 is not included in the I/O folder, 
this button still works), and it writes the following files, whose content will be 
described later: 

• probsF.txt 

• probsT.txt 

• blankets.txt 

• nits.txt 

Number of Probe Bits (for each PE step): This is the parameter a = 1, 2, 3, . . . 
for the operator V. 

Number of Phase Estimation (PE) Steps: This is the parameter c = 1, 2, 3, . . . 

for the operator V. 

McLximum Number of Grover Steps: Quibbs will stop iterating the AFGA if it 
reaches this number of iterations. 

Gamma Tolerance (degs): This is an angle given in degrees. Quibbs will stop 
iterating the AFGA if the absolute value of 7^- becomes smaller than this toler- 
ance. (7j is an angle in AFGA that tends to zero as the iteration index j tends 
to infinity. 7^ quantifies how close the AFGA is to reaching the target state). 

Delta Lambda (degs): This is the angle AA of AFGA, given in degrees. 

Once Quibbs has successfully read files Inl, In2 and In3, and once you have 
filled all the text boxes in the Inputs grouping, you can successfully press the Write 
Q. Circuit Files button. This will cause Quibbs to write the following output files 
within the I/O folder: 

(Outl) quibbs_log.txt 
(Out2) quibbs_eng.txt 
(Out3) quibbs_pic . txt 

The contents of these 3 output files will be described in detail in the next section. For 
this section, all you need to know is that: The quibbs_log.txt file records all the 
input and output parameters that you entered into the Control Panel, so you won't 
forget them. The quibbs_eng.txt file is an in"english" description of a quantum 
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circuit. And the quibbs_pic.txt file translates, line for line, the english description 
found in quibbs_eng.txt into a "pictorial" description. 

Normally, you want to press the Write Q. Circuit Files button without 
check-marking the Omit V Gates (diagnostic) check box. If you do check-mark 
it, you will still generate files 0utl,0ut2, Out3, except that those files will omit to 
mention all those gates that generate the operator V, at every place were it would 
normally appear. Viewing the circuit without its V^'s is useful for diagnostic and 
educational purposes, but such a circuit is of course useless for Gibbs sampling the 
Bayesian net being considered. 

The Control Panel displays the following output text boxes. (The Starting 
Gamma (degs) output text box and the Prob. of Starting State output text 
box are both filled as soon as a starting state is given in the inputs. The other output 
text boxes are filled when you press the Write Q. Circuit Files button.) 

Starting Gamma (degs): This is 70, the first 7j in AFGA, given in degrees. In 
the notations of Ref . [6] , and [8] , 



7o = acos(s' ■ t) = 2 acos(| I) (la) 
= 2 acos(/PO^) , (lb) 

where P{xo) is the Prob. of Starting State defined next. 

Prob. of Starting State: This is the probability P{xo) in Ref.[S], where P{) is the 
full probability distribution of the Bayesian net x being considered, and xq is 
the starting value for x. 

Number of Qubits: This is the total number of qubits used by the quantum circuit, 
equal to 2Nb + ac in the notation of Ref. [6] . 

Number of Elementary Operations: This is the number of elementary opera- 
tions in the output quantum circuit. If there are no LOOPs, this is the number 



of lines in the English File (see Sec. 4.2), which equals the number of lines in 



the Picture File (see Sec. 4.3). For a LOOP (assuming it is not nested inside 
a larger LOOP), the "LOOP k REPS:A^" and "NEXT k" hues are not counted, 
whereas the hues between "LOOP k REPS:A^" and "NEXT k" are counted 
times (because REPS:A^ indicates A^ repetitions of the loop body). Multiplexors 
expressed as a single line are counted as a single elementary operation (unless, 
of course, they are inside a LOOP, in which case they are counted as many 
times as the loop body is repeated). 
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3 Input Files 



As explained earlier, for Quibbs to generate quantum circuit files, it needs to first read 
3 input files: the Parents File called parents.txt, the States File called states.txt, 
and the Probabilities File called probs . txt. These 3 input files must be placed inside 
the I/O folder. Next we explain the contents of each of these 3 input files. 

3.1 Parents File 



© 3l parents.txt 




# A 




# B 




A 




C 




t C 









Figure 2: Parents file in the I/O folder "3nodes", for a Bayesian net with graph 

Fig|2] shows the Parents File as found in the folder "3nodes" which is included 
with the Quibbs distribution, for a Bayesian net with graph A B C. In this 
example, nodes A and C have no parents and node B has parents A and C. 

In general, a Parents File must obey the following rules: 

• Call focus nodes the node names immediately after a hash. Focus nodes in 
the States, Parents and Probabilities Files must all be in the same order. For 
example, in the "3nodes" case, that order is A, B, C. 

• For each focus node, give a hash, then the name of the focus node, then a list 
of parents of the focus node, separating all of these with whitespace. 

3.2 States File 

Fig. [3] shows the States File as found in the folder "3nodes" which is included with 
the Quibbs distribution, for a Bayesian net with graph A ^ B ^ C . In this example, 
node A has 3 states called al,a2 and a3, node B has 2 states called 61 and 62, and 
node C has 2 states called cl and c2. 

In general, a States File must obey the following rules: 

• Call focus nodes the node names immediately after a hash. Focus nodes in 
the States, Parents and Probabilities Files must all be in the same order. For 
example, in the "3nodes" case, that order is A, B, C. 
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MHIto _ 5tates;,.tJ<t 


# 


A 






al 








al 
















# 


B 






bl 








bl 








# 


c 






cl 








cl 






i 



Figure 3: States file in the I/O folder "Snodes", for a Bayesian net with graph A — > 
B^C 

• For each focus node, give a hash, then the name of the focus node, then a list of 
names of the states of the focus node, separating all of these with whitespace. 

3.3 Probabilities File 

Figjl] shows the Probabilities File as found in the folder "Snodes" which is included 
with the Quibbs distribution, for a Bayesian net with graph A B ^ C. In this 
example, PA(al) = 0.2, Pb|^ (7(6110!, cl) = 0.7, etc. 

In general, a Probabilities File must obey the following rules: 

• Call focus nodes the node names immediately after a hash. Focus nodes in 
the States, Parents and Probabilities Files must all be in the same order. For 
example, in the "Snodes" case, that order is A, B, C. 

• For each focus node, give a hash, then the name of the focus node, then the 
state of the focus node, then the states of each parent of the focus node, then the 
conditional probability of the focus node conditioned on its parents, separating 
all of these with whitespace. 

• The order in which the states of the parents of the focus node are listed must 
be identical to the order in which the parents of that focus node are listed in 
the Parents File. For example, in the "Snodes" case, the Parents File gives the 
parents of node B as A, C, in that order. Hence, in the Probabilities File, each 
conditional probability for focus node B is given after giving the states of nodes 
B, A, C, in that order. 

• A combination of node states may be omitted, in which case Quibbs will inter- 
pret that probability to be zero. For example, in "Snodes" case, if 
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© G JLI probs.txt 



# 


A 






Ql 








aZ 


0.3 






a3 
f^ 


0.5 
B 






bl 


Cll 


cl 


0.7 


bZ 


al 


cl 


0.3 


bl 


aZ 


cl 


0.4 


bZ 


aZ 


cl 


0.6 


bl 


a3 


cl 


0.1 


bZ 


a3 


cl 


0.9 


bl 


al 


cZ 


0.9 


bZ 


al 


cZ 


0.1 


bl 


aZ 


cZ 


0.5 


bZ 


qZ 


cZ 


0.5 


bl 


ai 


cl 


0.5 


bZ 


a3 


cZ 


0.5 


# 


C 






cl 


0.65 




cl 


0.35 





Figure 4: Probabilities File in the I/O folder "3nodes", for a Bayesian net with graph 
A-^ B^C 



PB|A,c(&l|a2,c2) = 0, 
you could omit a line of the form 

bl a2 c2 0.0 



for the focus node B. 

Note that Quibbs can help you to write a Probabilities File, by generating a 
template that you can change according to your needs. Such templates can be gener- 



ated by means of the Do a Pre-run button. See Sec j4.4| for a detailed explanation 
of this. 



4 Output Files 

As explained earlier, when you press the Write Q. Circuit Files button, Quibbs 
writes 3 output files within the I/O folder: a Log File called quibbs_log.txt, an 
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English File called quibbs_eng.txt, and a Picture File called quibbs_pic.txt. Next 
we explain the contents of each of these 3 output files. We also explain the contents 
of the various files generated when you press the Do a Pre-run button. 

4.1 Log File 



9 O 6 \=i quibbs_log,txt 

Inputs : 

I/O Folder = 3nodes/ 
Starting State : 
A; aZC01)l 
B; bZCl)l 
C: c 1(0)0 

Nymber Of Probe Bits (for each PE step) = 3 
Number Of Phase Estimation (PE) Steps = 3 
Maximum Number of G rover Steps = 10 
Gatrnia Tolerance Cdegs) =0.5 
Delta Lambda (degs) = 15^.00900006090003 

Outputs : 

3nodes/quibbs_eng.txt 
3nodes/quibbs_pic . txt 
3nodes/quibbs_log .txt 

Starting Gamma (degs) = 139,99603865314905 
Prob, of Starting State = 0,11699999999999999 
Wumber of Qubits = 17 
Number of Elem, Ops, = 8734 



Figure 5: Log File generated by Quibbs using input files from the "Snodes" I/O 
folder. 

FigjS] is an example a Log File. This example was generated by Quibbs using 
the input files from the "Snodes" I/O folder. A Log File records all the information 
found in the Control Panel. 

4.2 English File 

Fig|6] (respectively, Figj?]) is an example of an English File. This example was gen- 
erated by Quibbs, with the Omit V Gates feature OFF (respectively, ON), in the 
same run as the Log File of Figjsj and using the input files from the "Snodes" I/O 
folder. An English File completely specifies the output SEO. It does so "in English", 
thus its name. Each line represents one elementary operation, and time increases as 
we move downwards. 

In general, an English File obeys the following rules: 
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see 



CD 



1 


HADZ 


AT 


G 
















2 


HADZ 


AT 


7 
















M 


3 


HADZ 


AT 


3 


















4 


LOOP 


3 


REPSl 


1 














5 


PHAS 


180 




AT 


13 IF 


IZF IIF 


10 F 9F 


ST 








5 


MP_Y 


AT 


IZ 


IF 


15 C0 


14(1 


13C2 


ST 


BY 


-3^.76 




7 


MP_Y 


AT 


11 


IF 


14C0 


13(1 


IZCZ 


aT 


BY 


-33, Zl' 




8 


MP_Y 


AT 


10 


IF 


13 C0 


1Z(1 


11C2 


3T 


BY 


-30,85: 




9 


MP_Y 


AT 


9 


IF 


1ZC0 


11(1 


10CZ 


8T 


BY 


-4Z,79. 




10 


MP_Y 


AT 


16 


IF 


15 CZ 


14(1 


13C0 


ST 


BY 


39,75Z: 




1 1 


MP_Y 


AT 


15 


IF 


14CZ 


13(1 


1ZC0 


8T 


BY 


33,Z10I 




12 


MP_Y 


AT 


14 


IF 


13 CZ 


1Z(1 


11C0 


ST 


BY 


30.S63: 




13 


MP_Y 


AT 


13 


IF 


IZiZ 


11(1 


10C0 


8T 


BY 


4Z,794: 




14 


PHAS 


180,0 


AT 


IZ IF 


16F 15F 


14F 13 F 


8T 








15 


MP_Y 


AT 


13 


IF 


IZCZ 


11(1 


10C0 


ST 


BY 


-4Z,79. 




16 


MP_Y 


AT 


14 


IF 


13 CZ 


1Z(1 


11C0 


3T 


BY 


-30, S6: 




17 


MP_Y 


AT 


15 


IF 


14CZ 


13(1 


1ZC0 


8T 


BY 


-33,Zli 




18 


MP_Y 


AT 


16 


IF 


15 CZ 


14(1 


13C0 


8T 


BY 


-39,76; 


w 
T 


■ ^ ^ 





Figure 6: English File generated by Quibbs (with the Omit V Gates feature OFF) 
in the same run as the Log File of FigjSj and using the input files from the "Snodes" 
I/O folder. Bottom of file is not visible. Right hand side of file is not visible. 



• Time grows as we move down the file. 

• Each row corresponds to one elementary operation. Each row starts with 4 
letters that indicate the type of elementary operation. 

• For a one-bit operation acting on a "target bit" a, the target bit a is given after 
the word AT. 

• If the one-bit operation is controlled, then the controls are indicated after the 
word IF. T and F stand for true and false, respectively. aT stands for a control 
Pi (a) = n{a) at bit a. aF stands for a control Po{<y) = n{a) at bit a. 

• "LOOP k REPS : A^" and "NEXT k" mark the beginning and end of repetitions, 
k labels the loop, k also equals the line-count number in the English file (first 
line is 0) of the hue "LOOP k REPS : A^" . 

• SWAP a (3 stands for the swap(i.e., exchange) operator E{a,(3) that swaps bits 
a and (3. 

• PHAS 9'^^^'^ stands for a phase factor e*^''""" lio . 
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quibbs_eng.txt 



CD 



1 


P0PH 


150 . 00000000000003 


AT 


S 


IF 


7F 


6F 


5F 


4F 


3F 


ZF 


IF 


0F 


2 


P0PH 


160.7Z000357aZ5104 


AT 


9 


IF 


l&F 


15T 


14F 


13T 


12 F 


IIF 


10F 




3 


raPH 


150 , 00000000000003 


AT 


S 


IF 


7F 


6F 


5F 


4F 


3F 


ZF 


IF 


0F 


4 


P0PH 


124.139S49Z0956Z04 


AT 


9 


IF 


15F 


15T 


14 F 


13T 


12 F 


IIF 


10F 




5 


P0PH 


150 . 00000000000003 


AT 


S 


IF 


7F 


6F 


5F 


4F 


3F 


ZF 


IF 


0F 


5 


P0PH 


-5,069049526475199 


AT 


9 


IF 


16F 


15T 


14F 


13T 


12 F 


IIP 


10F 




1 


P0PH 


150 . 00000000000003 


AT 


S 


IF 


7F 


6F 


5F 


4F 


3F 


ZF 


IF 


0F 


8 


P0PH 


5,7339346S7546663 


AT 


9 


IF 


16F 


15T 


14 F 


13T 


12 F 


IIF 


10F 




9 


P0PH 


150 . 00000000000003 


AT 


S 


IF 


7F 


6F 


5F 


4F 


BP 


ZF 


IF 


0F 


10 


P0PH 


-3.897431743440155 


AT 


9 


IF 


15F 


15T 


14F 


13T 


12 F 


IIF 


10F 




1 1 


raPH 


150 , 00000000000003 


AT 


S 


IF 


7F 


6F 


5F 


4F 


3F 


ZF 


IF 


0F 


12 


P0PH 


4. 12 15 7 7035 110441 


AT 


9 


IF 


16F 


15T 


14 F 


13T 


12 F 


11 F 


10F 




13 


P0PH 


150 . 00000000000003 


AT 


S 


IF 


7F 


6F 


5F 


4F 


3F 


ZF 


IF 


0F 


14 


P0PH 


-Z,9S3S 843065765 84 


AT 


9 


IF 


16F 


15T 


14 F 


13T 


12 F 


IIP 


10F 




15 


P0PH 


150 . 00000000000003 


AT 


8 


IF 


7F 


6F 


5F 


4F 


3F 


2F 


IF 


0F 


16 


P0PH 


3,00160096716247 


AT 


9 


IF 


16F 


15T 


14 F 


13T 


12 F 


IIF 


10F 




17 


P0PH 


150 . 00000000000003 


AT 


S 


IF 


7F 


6F 


5F 


4F 


BP 


ZF 


IF 


0F 


18 


TOPH 


-Z.275Z69739356711S 


AT 


9 


IF 


16F 


15T 


14F 


13T 


12 F 


IIF 


10F 




19 


P0PH 


150 , 00000000000003 


AT 


S 


IF 


7F 


6F 


5F 


4F 


3F 


ZF 


IF 


0F 


20 


P0PH 


2,204S5S54S691914 


AT 


9 


IF 


16F 


15T 


14 F 


13T 


12 F 


11 F 


10F 




21 


P0PH 


150 . 00000000000003 


AT 


a 


IF 


7F 


6F 


5F 


4F 


3F 


ZF 


IF 


0F 


22 


P0PH 


-1,7ZSS9S096319023 


AT 


9 


IF 


16F 


15T 


14F 


13T 


12 F 


11 F 


10F 





Figure 7: English File generated by Quibbs (with the Omit V Gates feature ON) 
in the same run as the Log File of FigjSj and using the input files from the "Snodes" 
I/O folder. Bottom of file is visible. Right hand side of file is visible. 



• POPH Q'^''^' stands for the one-bit gate e*^"^''''" lio (note = n). PIPH Q'^''^' 
stands for the one-bit gate e^^^^ ""^iso (note P\ = n). Target bit follows the 
word AT. 

• SIGX, SIGY, SIGZ, HAD2 stand for the Pauli matrices ax-, <^y, <^z and the one-bit 
Hadamard matrix H, respectively. Target bit follows the word AT. 

• ROTX, ROTY, ROTZ, ROTN stand for one-bit rotations with rotation axes in the 
directions: x, y, z, and an arbitrary direction n, respectively. Rotation angles 
(in degrees) follow the words ROTX, ROTY, ROTZ, ROTN. Target bit follows the 
word AT. 

• MP_Y stands for a multiplexor which performs a one-bit rotation of a target bit 
about the y axis. Target bit follows the word AT. Rotation angles (in degrees) 
follow the word BY. Multiplexor controls are specified by a{k, where integer a 
is the bit position and integer k is the control's name. 
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Here is a list of examples showing how to translate the mathematical notation 
used in Ref. [Hj into the English File language: 



Mathematical language 


English File language 


J_jUUp licLilicU. O Willi Z( icpcLlLlUlio 


J_iUUr O riEjr O . Z 


INc-A-L ILcicLLlUli Ul lUUp licLllltjU. O 


IM I_iA i O 




ClUAP 1 n TF 9T 
owiir 1 \j ±r or zi 


„i42.7^n(3)n(2) 

±oU ^ ' ^ ' 


riiiio ^tz . f ±r or zi 


„i42.7^n(3)n(2) 


PDPH AO 7 AT TP OT 
ryjrn ^tZ . f iii O Xr Zi 


„i42.7^n(3)n(2) 


P1PH AO 7 AT TP OT 
rlrn ^tZ . 1 iii O Xr Zi 


n-^n'l"(3)"(2) 


CJTfY AT 1 TP '^P OT 

o±\jA iii 1 ±r or Zi 


o y 1 1 1 


cjTrv AT 1 TP '^P OT 
o±vjri iii 1 ±r or Zi 


n-^n ^"(3)n(2) 


CITny AT 1 TP '^P OT 
ij±\3£j iii 1 xr or Zi 


i^(l)n(3)n(2) 


HAD2 AT 1 IF 3F 2T 


(e*lfo23-7<Tx{l))n(3)n(2) 


ROTX 23.7 AT 1 IF 3F 2T 


(e*lfo23-7<Ty(l))n(3)n(2) 


ROTY 23.7 AT 1 IF 3F 2T 


(e*iio23-7<^z(i))"(3)"(2) 


ROTZ 23.7 AT 1 IF 3F 2T 


(^gj^[30ax(l)+40<7y(l)+lW(l)]^|n{3)n(2) 


DdTM Qn n /in n -i 1 n at -i tt? qt? ot 
rtUllN OU . U ^tU . U ll.U Ai 1 Ir or Zi 


j^g'Ebl.bO 

where < 


' 000 = 30.0(^) 
^^01 = 10.5(^) 
9^o = ll.O(^) 

^ 011 = 83.1(jfo) 


MP_Y AT 3 IF 2(1 1(0 OT BY 30.0 10.5 11.0 83.1 



4.3 ASCII Picture File 

Figjs] (respectively, Fig|9]) is an example of a Picture File. This example was generated 
by Quibbs, with the Omit V Gates feature OFF (respectively, ON), in the same 
run as the Log File of Figjsj and using the input files from the "Snodes" I/O folder. 
A Picture File partially specifies the output SEO. It gives an ASCII picture of the 
quantum circuit. Each line represents one elementary operation, and time increases 
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eoe 


quibbs_pic.txt 


CD' 




1 I 


1 1 1 1 1 1 1 1 H 1 1 1 




m 
1 


2 


1 1 


1 1 1 1 1 1 1 H 1 1 1 1 


1 1 1 


V 


3 


1 1 


1 1 1 1 1 1 H 1 1 1 1 1 


1 1 1 




4 


LOOP 


3 REPS; L 






5 


1 1 


1 ph — — — — — 1 1 1 1 1 


1 1 1 




6 


1 (0- 


-CI — Q.~-^'^-—v — + — + — 'S 1 1 1 1 1 


1 1 1 




1 


1 1 


^0__C1--(Z--Ry--+ +---@ 1 1 1 1 1 


1 1 1 




8 


1 1 


1 -f 1 - -f 7 - -Rv- 1 1 1 1 1 


1 1 1 




Q 


1 1 


1 1 C0--Ci--C^--Ry--'^ 1 1 1 1 1 






IS 


Ry--CZ- 


-C1--C0 — 1 1— ..+.-^+^ — % 1 1 1 1 1 






1 1 


1 Ry- 


-C2—C1—C0- -+---+— 1 1 1 1 1 


1 1 1 




12 


1 1 


Ky — — s^l--\_^a 1 + ^ 1 1 1 1 1 






13 


1 1 


1 f? n f'Ck J. S 1 1 III 
1 Ky — \£. — \l. — \_\f 1 1 1 1 1 1 






14 


0---0— 


-0-— 0---Ph--+---+— -+— @ 1 1 1 1 1 






15 


1 1 


1 Ry--(Z--Cl--C0- -+-—'§ 1 1 1 1 1 






16 


1 1 


Ry--CZ--Cl--C0--+"-+— ^ 1 1 1 1 1 


1 1 1 




17 


1 Ry- 


-CZ— C1--C0--+---+— @ 1 1 1 1 1 




A. 
T 


o 











Figure 8: Picture File generated by Quibbs (with the Omit V Gates feature OFF) 
in the same run as the Log File of Figjsj and using input files from the "Snodes" I/O 
folder. Bottom of file is not visible. 



as we move downwards. There is a one-to-one onto correspondence between the rows 
of the English and Picture Files. 

In general, a Picture File obeys the following rules: 

• Time grows as we move down the file. 

• Each row corresponds to one elementary operation. Columns 1, 5, 9, 13, . . . rep- 
resent qubits (or, qubit positions). We define the rightmost qubit as 0. The 
qubit immediately to the left of the rightmost qubit is 1, etc. For a one-bit 
operator acting on a "target bit" a, one places a symbol of the operator at bit 
position a. 

• I represents a "qubit wordline" connecting the same qubit at two consecutive 
times. 

• -represents a wire connecting different qubits at the same time. 

• + represents both I and -. 

• If the one-bit operation is controlled, then the controls are indicated as follows. 
at bit position a stands for a control n{a) = Pi (a). at bit position a stands 
for a control n{a) = Po{a)- 
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©e© 








ji^ quibbs_pic.txt 












CD 




1 


1 1 1 


1 


1 


1 1 


0P--O---O- 


._o 


Q__ 


-0- 


--0- 


--0 






2 






--0- 




-0— -OP 


1 1 1 


1 


1 1 


1 


1 


1 






3 


1 


1 1 1 


1 


1 


1 1 


0P--O---O- 


- -0-- 


Q__ 


-0- 


--0- 


- -0 






4 


0--- 


a Q @_ 


Q_ 


--0-- 


-0— -OP 


1 1 1 


1 


1 1 


1 


1 


1 






5 


1 


1 1 1 


1 


1 


1 1 


OP--0 -0- 


-_Q 


-Q Q__ 


-0- 


--0- 


--0 






5 


0--- 


a Q ig. 


- -0- 


--0-- 


-0— -OP 


1 1 1 


1 


1 1 


1 


1 


1 




1 


1 


1 


1 1 1 


1 


1 


1 1 


0P--O---O- 


._Q 


-Q, Q__ 


-0- 


--0- 


--Q 






8 


0--- 


a Q ^_ 


--0- 


--0-- 


-0— -OP 


1 1 1 


1 


1 1 


1 


1 


1 






9 


1 


1 1 1 


1 


1 


1 1 


0P--O---O- 


--0 





-0- 


--0- 


--0 






10 


0--- 


a Q @_ 


--0- 


--0-- 


-0— -OP 


1 1 1 


1 


1 1 


1 


1 


1 






1 1 


1 


1 1 1 


1 


1 


1 1 


0P--O---O- 


._Q 


Q__ 


-0- 


--0- 


--0 






12 


0--- 


a Q @_ 


--0- 


--0-- 


-0— -OP 


1 1 1 


1 


1 1 


1 


1 


1 






13 


1 


1 1 1 


1 


1 


1 1 


0P--O---O- 


- -0-- 


-Q Q__ 


-0- 


--0- 


- -0 






14 


0--- 


a Q g_ 


- -0- 


--0-- 


-0— -OP 


1 1 1 


1 


1 1 


1 


1 


1 






15 


1 


1 1 1 
1 1 1 


1 

1 


1 

1 


1 1 


0P--O---O- 


- -fl-- 


_n n__ 












1 D 


u 


e>---0---@'- 


— 0- 


--0-- 


n. no 

-\j Ur 


1 1 1 
1 1 1 


1 


1 1 


1 


1 


1 






\ 7 


1 

1 


1 1 1 


1 


1 


1 1 
1 1 


np- -O ^O- 


--0— 


.0—0-- 


-0- 


--0- 


--0 






18 


0--- 


e— -0— @- 


— 0- 


--0-- 


-0— -OP 


1 1 1 


1 


1 1 


1 


1 


1 






19 


1 


1 1 1 


1 


1 


1 1 


0P--O---O- 


--0— 


-0— 0-- 


-0- 


--0- 


--0 






20 


0--- 


e---o---@'- 


— 0- 


-0- 


-0— -OP 


1 1 1 


1 


1 1 


1 


1 


1 




J 


21 


1 


1 1 1 


1 


1 


1 1 


0P--O---O- 


--0— 


.0—0-- 


-0- 


--0- 


--0 






22 


0--- 


g— -0— g- 


--0- 


--0-- 


-0 — -OP 


1 1 1 


1 


1 1 


1 


1 


1 




X 

T 




























/a 



Figure 9: Picture File generated by Quibbs (with the Omit V Gates feature ON) 
in the same run as the Log File of Figjsj and using input files from the "3nodes" 1/0 
folder. Bottom of file is visible. 



• "LOOP k REPS : A^" and "NEXT k" mark the beginning and end of N repetitions, 
k labels the loop, k also equals the line-count number in the Picture File (first 
line is 0) of the line "LOOP k REPS : A^" . 

• The swap(i.e., exchange) operator £'(«, /3) is represented by putting arrow heads 
< and > at bit positions a and /3. 

• A phase factor e*^ for ^ G M is represented by placing Ph at any bit position 
which does not already hold a control. 

• The one-bit gate e*^o(")^ (note P^ioL) = n{a)) for ^ G M is represented by putting 
OP at bit position a. 

• The one-bit gate e*^^^"^^ (note Pi{a) = n{a)) for 6' G M is represented by putting 
OP at bit position a. 

• One-bit operations ax (a), cry (a), cr^(a) and H{a) are represented by placing 
the letters X,Y,Z, H, respectively, at bit position a. 
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• One-bit rotations acting on bit a, in the x, y, n directions, are represented by 
placing Rx,Ry,Rz, R, respectively, at bit position a. 

• A multiplexor that rotates a bit r about the y axis is represented by placing 
Ry at bit position r. A multiplexor control at bit position a and named by the 
integer k is represented by placing {k at bit position a. 

Here is a list of examples showing how to translate the mathematical notation 
used in Ref. [9] into the Picture File language: 



Mathematical language 


Picture File language 


Loop named 5 with 2 repetitions 


LOOP 5 REPS: 2 


Next iteration of loop named 5 


NEXT 5 


£;(l,0f(3M2) 


— @ — < — > 


gi42.7^n(3)n{2) 


0— @— +— Ph 


gi42.7^n(3)n(2) 


OP— @ 1 




gi42.7^n(3)n{2) 


@P— @ 1 




^^^l)n(3)n(2) 


0— @— X 




ay(l)"(3)"(2) 


0— @— Y 




^^(l)n(3)n(2) 


0— @— Z 




if(l)"(3)n(2) 


0— @— H 




(^gj^23.7<7x{l))n(3)n{2) 


@ Rx 




^gi^23.7<7y(l)^)n(3)n(2) 


0— @ Ry 




(e*lio23-7'Tz(l))n(3)n(2) 


@ Rz 




(-gi^[30<7x(l)+40ay(l)+llaz(l)]^|n{3)n(2) 


0— @— R 




[e'Ebi.bo Sbo'^i'(3)P6^6o{2,l)jn(0) 


1 Ry— (1— (0— @ 



' 000 = 30.0(^) 



^^oi = 10.5(^) 
^io = 11.0(^) 
^ 9n = 83.1(^) 
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4.4 Output Files From a Pre-run 



When you press the Do a Pre-run button, Quibbs writes 4 output files within the 
I/O folder: two Uniform Probabilities Files called probsF.txt and probsT.txt, a 
Blankets File called blankets .txt, and a Nits File called nits .txt. Next we explain 
the contents of each of these 4 output files. 

4.4.1 Uniform Probabilities Files 

IjpPBWP [g probsF.txt 

# A 

a5 0.3333333333333333 



bl 


al 


cl 


0.5 


bl 


al 


cl 


0.5 


bl 


aZ 


cl 


0,5 


bl 


aZ 


cl 


0.5 


bl 


a3 


cl 


0.5 


bl 


a3 


cl 


0.5 


bl 


al 


tl 


0.5 


bl 


al 


cl 


0.5 


bl 


aZ 


cZ 


0.5 


bl 


aZ 


cZ 


0,5 


bl 


a3 


cZ 


0.5 


bl 


a3 


cZ 


0.5 


# 


C 






cl 


0,5 






a 


0.5 







Figure 10: Uniform Probabilities File (of type probsF.txt) generated using input 
files from the "3nodes" I/O folder. 



FigsjTO] and 11 are both examples of a Uniform Probabilities File. Both files 
were generated by Quibbs using the input files from the "Snodes" I/O folder. A 
Uniform Probabilities File is simply a Probability File, as defined in Sec, 3.3, but of a 
specific kind that assigns uniform values to all conditional probabilities of the Bayesian 
net specified by the Parents File and States File in the I/O Folder. Generating 
a Uniform Probabilities File (either probsF.txt or probsT.txt) does not require 
a probs . txt file. One can use a Uniform Probabilities File as a template for a 
probs .txt file. Just cut and paste the contents of a Uniform Probabilities File into a 
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© G S probsT.txt 



# 


A 






al 


0,3333333333333333 


aZ 


0,3333333333333333 




0.3333333333333333 


# 


B 






bl 


Ql 


cl 


0.5 


bl 


aZ 


cl 


0.5 


bl 


0.3 


cl 


0.5 


bl 


Ql 


cZ 


0.5 


bl 


qZ 


cZ 


0.5 


bl 


0.3 


cZ 


0.5 


bZ 


Ql 


cl 


0.5 


bl 


qZ 


cl 


0.5 


bl 


03 


cl 


0.5 


bZ 


ol 


cZ 


0.5 


bZ 


aZ 


cZ 


0.5 


bl 


03 


cZ 


0.5 


# 


C 






cl 


0.5 






c2 


0,5 







Figure 11: Uniform Probabilities File (of type probsT.txt) generated using input 
files from the "3nodes" I/O folder. 



new file called probs . txt and modify its probabilities according to your needs. Note 
that the only difference between a probsF . txt and a probsT . txt file is that the first 
(respectively, second) of these varies the states of a focus node before (respectively, 
after) varying the states of its parents. 

4.4.2 Blankets File 



Figjl2]is an example of a Blankets File. This example was generated by Quibbs using 
the input files from the "3nodes" I/O folder. 

The Markov blanket MB{i) for a node of the classical Bayesian network x 
is defined so that (see section entitled "Notation and Preliminaries" in Ref.|lU]) 

P{Xi\x{iY) = P{Xi\xMB{i)) ■ (2) 

It can be shown that the Markov blanket of a focus node equals the union of: 

• the parents of the focus node 

• the children of the focus node 



18 



© G blanket5.txt 



# 


A 




a 






c 






# 


B 




A 






c 






# 


C 




A 













Figure 12: Blankets File generated using input files from the "3nodes" I/O folder. 



the parents of each children of the focus node (but excluding the focus node 
itself) 

A Blankets File gives the Markov blanket for each node of a Bayesian net. In 



the example of Fig, 12, node A has Markov blanket {-B, C}, etc. 



In general, a Blankets File obeys the following rules: 

• Call focus nodes the node names immediately after a hash. 

• For each focus node, Quibbs writes a hash, then the name of the focus node, 
then a list of the nodes which form the Markov blanket of the focus node. 

4.4.3 Nits File 



Fig, 13 is an example of a Nits File. This example was generated by Quibbs using the 
input files from the "Snodes" I/O folder. 

The word "nit" is a contraction of the words "node" and "qubit". Quibbs 
assigns to each node (of the Bayesian net being consider) its own private set of nits. 
We explained in Sec|2]how Quibbs assigns a decimal and a binary name to each state 
of a node. The binary name of a state gives the states of the nits. For example, 
suppose node A has 3 states: al(00)0, al(01)l and al(10)2. Then node A is assigned 
two private nits, call them nit^ and n%t\. When node A is in state (61, 60), where 61 
and 60 are either or 1, then nit\ is in state 61 and nitO is in state 60. 

Actually, Quibbs doesn't give nits an "english" name like nitO and n%t\. It 



just calls them by integers. Fig,13 informs us that the "3nodes" Bayesian net has 
4 nits called 0,1,2,3. Nits and 1 are both owned by node A (which has 3 states 
al, a2, a3). Nit 2 is owned by node B (which has 2 states 61, 62). Nit 3 is owned by 
node C (which has 2 states cl, c2). 

The original Bayesian net with the original nodes implies a new, finer Bayesian 
net whose nodes are the nits themselves. Just like one can define a Markov blanket 
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© o e 


nits.txt 






# 










ownerMode: A 










nit blanket; 


(1 in A)CZ in B)(3 


in 


C) 




# 1 










cwvnerMode: A 










nit blanket; 


(0 in A)CZ in B)(3 


in 


C) 




# Z 










ownerMode; B 










nit blanket; 


(0 in A)C1 in AX3 


in 


C) 




# 3 










ownerMode; C 








1 


nit blanket ; 


(0 in AXl in AX^ 


in 














A 



Figure 13: Nits File generated using input files from the "3nodes" I/O folder. 



for each node of the original Bayesian net, one can define a Markov blanket ( equal to 



a particular set of nits) for each nit. Fig, 13 informs us that for the "3nodes" Bayesian 
net, nit has a nit blanket {1, 2, 3}, etc. 

In general, a Nits File obeys the following rules: 

• Call focus nit the number ( a sort of nit name) immediately after a hash. 

• For each focus nit, Quibbs writes the words "owner node" followed by the name 
of the owner node of the focus nit. 

• For each focus nit, Quibbs writes the words "blanket nit" followed by the Markov 
blanket of nits for the focus nit. 

Why do we care about node blankets and nit blankets at all? Quibbs uses 
a method, discussed in Ref. [;6], of representing Szegedy operators using quantum 
multiplexors. Quibbs uses nit blankets to simplify its Szegedy representations by 
eliminating certain unnecessary controls in their multiplexors. 
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